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A weakly coupled scalar field <E> with a simple exponential potential V = Mp exp(— \&/Mp) 
where Mp is the reduced Planck mass, and A > 2, has an attractor solution in a radiation or matter 
dominated universe in which it mimics the scaling of the dominant component, contributing a fixed 
fraction (determined by A) to the energy density. Such fields arise generically in particle physics 
theories involving compactified dimensions, with values of A which give a cosmologically relevant 
£1$. For natural initial conditions on the scalar field in the early universe the attractor solution 
is established long before the epoch of structure formation, and in contrast to the solutions used 
in other scalar field cosmologies, it is one which does not involve an energy scale for the scalar 
field characteristic of late times . We study in some detail the evolution of matter and radiation 
perturbations in a standard inflation-motivated Q = 1 dark-matter dominated cosmology with 
this extra field. Using a full Einstein-Boltzmann calculation we compare observable quantities with 
current data. We find that, for fi^, ~ 0.08 — 0.12, these models are consistent with large angle cosmic 
microwave background anisotropies as detected by COBE, the linear mass variance as compiled from 
galaxy surveys, big bang nucleosynthesis, the abundance of rich clusters and constraints from the 
Lyman-a systems at high redshift. Given the simplicity of the model, its theoretical motivation and 
its success in matching observations, we argue that it should be taken on a par with other currently 
viable models of structure formation 

PACS Numbers : 98.80.Cq, 98.70.Vc, 98.80.Hw 



I. INTRODUCTION 

The past twenty years have seen a tremendous revolution in how we study the origin and evolution of our universe. 
On the one hand developments in theoretical particle physics have lead to a proliferation of ideas on how high energy 
physics might have an observable effect on the large scale structure of our universe. On the other hand the increasing 
quality of astrophysical data has led to firm constraints on what physics is allowed in the early universe. Probably 
the most impressive example of such an interplay is how the COBE detection jjj has affected the most popular and 
theoretically explored theory of structure formation, the standard cold dark matter model (SCDM). 

The SCDM model brings together the idea of inflation ^j-Q and the picture of large scale gravitational collapse 
||. A period of superluminal expansion of the universe would have led to the amplification of subhorizon vacuum 
fluctuations to superhorizon scales. The net result would be a set of scale invariant, Gaussian perturbations which 
would evolve, through gravitational collapse, into the structures we see today. The relic radiation bears an imprint 
of these fluctuations from when the universe recombined. The distribution of galaxies and clusters of galaxies should 
reflect these fluctuations today. It has been found however that the SCDM model cannot successfully accommodate 
the data observed on all scales. Matching its predictions to COBE measurements (on large scales) of the microwave 
background, one finds that the amplitude of fluctuations on 8h _1 Mpc scales (where h is the Hubble constant today 
in units of lOOkm/s/Mpc ) do not match those observed jfj). To be more specific, if we define the amplitude of mass 
fluctuations in a sphere of radius 8h _1 Mpc, as, then COBE normalized SCDM predicts erg = 1.2 while the measured 
value (through cluster abundances) is as — 0.6 ± 0.1, a discrepancy by a factor of 2. Further not only the amplitude 
but also the scale dependence of the SCDM model differs from the one measured 0, and there are also a number of 
problems with the non-linear evolution of baryons and velocities on small scales. 

These failings of SCDM have led to attempts to modify it, while keeping its basic features intact. The latter (a simple 
choice of background cosmological parameters with 17 = 1 and quantum generation of fluctuations) are associated with 
its theoretical motivation from inflation. The most prominent candidate theories of structure formation of this type are 
now a universe with a cosmological constant (ACDM) ]8| and a universe with a fifth of its dark matter component in 
two families of massive neutrinos (MDM) Q. Models of structure formation in an open universe (OCDM), for which 
there is considerable evidence, have also been extensively studied jl(|. All these models, once COBE normalized, 
predict the approximately correct distribution of mass fluctuations. Over the past few years other flat models have 
been constructed which, like MDM, are more related to our understanding of fundamental particle physics. This is 
the case of decaying cosmological constant models [jl2],[l3| and decaying massive particles |]|l4|] • Unfortunately, unlike 
the SCDM scenario, all these models involve a tuning of parameters which is unnatural from the point of view of 
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particle physics, simply because one is using super-GeV physics to explain sub-eV observations. Just as is the case of 
ACDM (which involves tuning the cosmological constant to be relevant only at present epochs) this tuning does not 
provide a reason to discard these models, but is a very unattractive feature of them. 

Both the cosmology of weakly coupled scalar fields and their theoretical motivation have been much studied since 
the advent of the idea of inflation. In most of the theoretical particle physics extensions of the standard model (e.g. 
supersymmctry and supergravity) weakly coupled scalar fields of various kinds are ubiquitous. Unfortunately there 
is as yet no experimental evidence for the existence of any fundamental scalar field at all, in particular the Higgs 
particle of the standard model remains to date undetected. However the theoretical motivation for such scalar fields 
is sufficiently compelling that it is certainly worthwhile considering what consequences their existence might have for 
cosmology beyond the confines of inflation. An example of this is given by the main alternative of structure formation 
- defect theories - which usually involve the existence of some scalar field e.g. the 'texture' theory jl5| relies on the 
existence of a scalar field invariant under a global non-abelian symmetry broken at the GUT scale, with associated 
goldstone particles which are unobservably weakly coupled to ordinary matter. 

The role that a weakly coupled scalar field might play in late time cosmology if it contributes a component to the 
homogeneous background energy density of the universe has also been investigated. In ||l6| the authors considered 
in general terms the idea that a significant contribution to the energy density from a homogeneous scalar field could 
have an important effect on structure formation, and applied the idea to a baryonic universe with a initial spectrum of 
perturbations with power law behaviour. After some general analysis of the homogeneous dynamics, the model singled 
out for a detailed treatment was that of a scalar field in a negative power law potential which comes to dominate at 
late times, producing behaviour very similar to that associated with a cosmological constant . Another kind of model 
which has been developed in detail, for the case of a cold dark matter dominated cosmology, in flj,[]l| involves a 
cosine potential with a combination of Planck and sub eV physics. This provides a specific realisation of a 'decaying 
cosmological constant', in which the field initially behaves like a cosmological constant and then rolls so that it scales 
asymptotically like matter. All these models have an energy density in the scalar field which comes to dominate at 
late times, as have two more recent very detailed studies |l7],[l8|]. In fTiJ the case of a scalar field which scales slower 
than matter is described more generically in terms of its equation of state, and in |isj | the specific cases of late time 
dominance realized in both cosine and exponential potentials. 

In this paper we present in detail the results which have been reported in The model we study is SCDM with 
the addition of a scalar field $ with a simple exponential potential V — Mp exp(— A$/Mp) where Mp = (8ttG)~z = 
2.4 x 10 18 GeV is the reduced Planck mass, and A > 2. In this specific case there is a very special and interesting 
solution for matter and radiation coupled through gravity to this field, in which the scalar field energy follows that of 
the dominant component, contributing a fixed fraction of the total energy density determined by A. The existence of 



this homogeneous solution was shown in |pL 6|J2 0|J2 If] . Because the scalar field can contribute at most a small fraction of 
the total energy density at nucleosynthesis, its potential interest in the context of the problem of structure formation 
has been overlooked. This constraint suggests prima facie that the field in such an attractor can have little effect on 
cosmology at late time. As we shall see, this is incorrect, for the simple reason that a small contribution acting over 
a long time can have as big an effect as a large contribution entering only at late times. The particular merit of the 
model is related to the fact that the cosmological solution is an attractor for the scalar field: Because of this, there is 
no tuning of the type involved in all other proposed modifications of SCDM. The only parameter additional to SCDM 
is A, and the value (~ 5 — 6) which gives a best fit to structure formation, is of the order naturally expected in the 
particle physics models in which the potential arises. As a cosmology it resembles MDM much more than any of the 
scalar field models which have been studied in the literature. 

Some comment is perhaps necessary at the outset regarding the assumption that the universe is flat and dominated 
by a component scaling as matter, as there is mounting observational evidence that this is not the case. The most 
advertised is the age problem, that is, the fact that we see objects which are older than the age of the universe if we 
assume a flat matter dominated universe with the currently observed Hubble constant of Ho = 65 ± 10kms _1 Mpc _1 . 
Several of the modifications we have discussed (in particular the cosmological constant model and some of the scalar 
field models) avoid this problem since the dominant component at the present epoch is not matter. The recent 
measurements of the Hipparcos satellite |22| seem to indicate that the age of these objects have been overestimated 
by around 10% and a re-analysis of the uncertainties in the estimates seem to indicate that a flat matter dominated 
universe with H$ < 66kms _1 Mpc _1 is compatible with the current age estimates. Another argument for a low density 
universe comes from the analysis of large scale flows. In |23 it has been argued that (3 = fi°- 6 /& < 1 from the analysis 
of the MARK III data and comparison with the IRAS surveys. However there is concern with the self consistency of 
the velocity data and another group |£4| indicates that the the flows are consistent with a universe with £1 = 0.4 — 1. 
Small scale observations also seem to indicate that there are problems with a high density universe. In particular the 
baryon fraction in clusters is difficult to reconcile with the BBN limits unless one considers a low density universe [ p6[ . 
Again there are uncertainties in such an analysis; they rely on elaborate numerical simulations which are at the limit 
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of current computational power. It is conceivable that some physics (such as cooling of the cluster medium) is being 
overlooked. The same can be said for the cold velocity dispersion which is measured on few Mpc scales ||. Although 
all these observations together begin to make a strong case for an open universe, the uncertainties and inconsistencies 
are sufficiently large for us to still consider a high density, matter dominated universe. It may be that in the next few 
years the evidence is sufficiently compelling to rule out these models, but this is definitely not the case yet. 

The structure of the paper is as follows. In section || we first discuss the homogeneous modes of scalar fields 
in a universe with matter and radiation in general terms, and explain how the case of a simple exponential with 
its attractor solutions for A > 2 is a special one. We discuss briefly the possible origin of exponential potentials 
in fundamental theories. We then discuss the initial conditions on this scalar field in the early universe, and when 
these lead the attractor to be established. In typical inflationary theories we argue that the attractor describes the 
homogeneous cosmology well before nucleosynthesis, while in an alternative theory of reheating which can be realised 
with the same exponential field the scalar field may still be negligible at nucleosynthesis for natural parameter values. 
In section [II we analyse the evolution of perturbations in our scenario. We describe the complete set of equations 
which govern their evolution and analyse the asymptotic behaviour in the interesting regimes. The similarities with 
evolution of density fluctuations in an MDM universe lead us to pursue the comparison in some detail and we come to 
the conclusion that scalar field is more effective at suppressing perturbations than hot dark matter. We also consider 
the effect on the cosmic microwave background (CMB) and deconstruct the different effects it has on the angular 
power spectrum of the CMB. In section IV we compare the predictions of linear perturbation theory (using the full 
Boltzmann-Einstein equations) with some observations. In particular we compare it to the COBE data and then use 
this comparison to normalize the theory and compare to the Peacock and Dodds data. As a byproduct we derive a 
fit to 08 for our models and we quantify the amount of structure on small scales and high redshift, comparing with 
constraints derived from Lyman-a systems. In ^ we summarise our findings and draw some conclusions about future 
prospects for our model and other issues related to it which might be further investigated. 



II. THE 'SELF-TUNING' SCALAR FIELD WITH AN EXPONENTIAL POTENTIAL 



In this section we describe in detail the homogeneous attractor solutions which specify the zero-order cosmology, 
about which we treat perturbations fully in the subsequent part of the paper. We explain the very special feature of 
this model which contrasts it with other scalar field cosmologies which have been treated in the literature: No special 
energy scale characteristic of late time cosmology need be invoked to produce the required solution. We will show that 
for a very wide and natural range of initial conditions on the scalar field in the very early universe, the attractor will 
be attained as early as assumed in the rest of the paper. 



A. Scaling of the Energy Density in a Scalar Field 



We begin with a general discussion of homogeneous FRW cosmology in which there is, in addition to the usual 
matter and radiation content, a contribution to the energy momentum coming from a scalar field <f> with a potential 
V(&). In the rest of the paper we use the notation $ = 4>+<p, where ip denotes the perturbation about the homogeneous 
solution <f>. Working for the moment in comoving coordinates in which the metric is ds 2 — —dt 2 + a 2 8ijdx l dx'^ (where 
a is the scale factor), the contribution to the energy- momentum tensor from the scalar field is 
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where 



The equations of motion for the scalar field are then 

<f> + 3H4> + V'(<f>) = 
H 2 = 



U 2 + V{cj>) 



1 ■ 



-J 2 - V(</>) 



nHp n = 



(1) 

(2) 

(3) 
(4) 



where p n is the energy density in radiation (n = 4) or non-relativistic matter (n = 3), H = — is the Hubble expansion 
rate of the universe, dots are derivatives with respect to time, primes derivatives with respect to the field <fi, and 
Alp = (8ttG)-2 = 2.4 x 10 18 GeV is the reduced Planck mass. The scalar field is assumed to be coupled to ordinary 
matter only through gravity. Multiplying (0) by (f> and integrating, one obtains 



3 



(a) = p(a a ) exp ( — f 6(1 — £(a))— ) where 



+ c = ^. (5) 



It follows from this that, given the range of possible values < £ < 1 (assuming V((p) is positive), the energy density 
of a scalar field has the range of scaling behaviours 

P4> (xl/a m 0<m<6 (6) 

How the energy in a homogeneous mode of a scalar field scales is thus determined by the ratio of the potential to 
the kinetic energy. Alternatively one can phrase the statement in terms of the equation of state obeyed by the mode: 
From (|l|) we have £ = i(l — to) where = wp^, and m = 3(1 + to) (for constant w) in (||). 

These statements are true independent of any specific assumption about H i.e. about what dominates the energy 
density of the universe. When the potential energy dominates over kinetic energy, we have £ — > 1 and therefore 
Prj, rs constant i.e. an energy density which behaves like a cosmological constant (and w — —1); in the opposite limit 
of £ — > i.e. a kinetic energy dominated mode, we have an energy density with p^ oc 1/a 6 i.e. red-shifting faster than 
radiation or matter (and to = +1). Inflation occurs when the former type of mode also dominates the energy density; 
the opposite limit, when the universe is dominated by the kinetic energy of a scalar field (which, following Q], we 
refer to as kination) gives a sub- luminal expansion with a oat's . 

We now consider more specifically what sorts of potentials give rise to these different types of scaling. One simple 
case to analyse is that in which a field rolls about the minimum of a potential. In such an oscillatory phase the 
approximate scaling of the energy density can be extracted from (||) by replacing £ by its average value over an 
oscillation. For a potential which is power law about its minimum with V(<p) oc cj> n the result | p7| is that 

m=^- 2 (7) 

reproducing the well known result that a coherent mode oscillating in a quadratic potential gives the same scaling as 
matter, and a 4 potential that of radiation. For n > 4 one obtains modes scaling faster than radiation. 

Again this statement does not depend on what component dominates the energy density, and the same scaling 
applies to the mode irrespective of whether the universe has p n = 0, or is matter or radiation dominated. The case 
of a field rolling down a potential (before it reaches its minimum, or if it has no minimum) is quite different. The 
equation of motion (j^) is just a damped roll with the energy content determining the damping through (j^). The 
scaling obtained for a given potential depends on what components are present, because what determines the scaling 
is the balance between the increase in kinetic energy relative to potential energy as the field rolls down the potential, 
and the decrease of the same quantity due to the damping. The criterion for a particular scaling is therefore a 
requirement of the 'steepness' of the potential. That the simple exponential potential V((f>) = V e^ x ^/ Mp provides 
the appropriate yard-stick in the case of scalar field domination is indicated by the existence of a family of solutions 
for this potential to @-@ with p n = @||>]|2l]] 

m=to+ 2 -^MtM P ) *o = ^H 2M ?f_ x2) ) e = l-y P<»oc± fl «^ 2 (8) 

for A < V6, and <j) can always be chosen to be zero by redefining the origin of (j>, in which case V = — l)Mp. 

These solutions, which are attractorsQ were written down ^8|,^9| in the context of power-law inflation p0[ |, associated 
with the superluminal growth of the scale-factor for A < v2. For A > 6 there is not a single attractor with a finite 
value of £, but every solution has £ — ► asymptotically and p oc 1/a 6 . 

As A increases from zero we obtain the entire range of possible scaling behaviours @ for the energy in a scalar field. 
By comparison with this potential we can infer how scalar modes will scale: The 'slow-roll' conditions (\MpV' /V\ << 
V6, \MpV" /V\ « 3) for inflation, for example, could be stated as the requirement that the first two derivatives of 
a potential be smaller than those of an exponential with A = An analagous 'fast-roll' condition for 'kination' 
can clearly be provided by comparison with an exponential with A = \/6 e.g. a potential ~ e - ^ will have such 
modes for sufficiently large <j>. 

These statements apply only to the case of scalar field dominance since we took p n = 0. What interests us in the 
present context is how the scalar energy behaves in the presence of matter and radiation i.e. with p n =/= 0. There are 



1 It is simple to verify directly from (g) and (^) that homogeneous perturbations about the solution decay as t 1 and t 1 
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FIG. 1. In the left panel we plot the evolution of the energy density in the scalar field (p^) and in a component of radia- 
tion-matter as a function of scale factor for a situation in which the scalar field (with A = 4) initially dominates, then undergoes 
a transient and finally locks on to the scaling solution. In the right panel we plot the evolution of the fractional density in the 
scalar field. 



two quite distinct cases which can be immediately distinguished according to their behaviour when p n = 0: Those 
potentials in which the energy density scales slower than 1/a™ (i.e. A < y/n), and those in which it scales faster 
(A > y/n). When p n ^ these two types of potential will show very different behaviour for the following reason: 
Adding a component increases the damping term in (|J), and it follows that the scaling with a of the energy density 
in the scalar field is always slower i.e. p^ oc l/a x ~ s with A 2 > S > 0. For A < y/n the scalar energy will still red-shift 
slower than the other component and it will always come to dominate asymptotically, approaching the attractor (^). 
For A > y/n, however, the scalar field energy cannot always scale as in the case p n = 0. By doing so it would 
become arbitrarily sub-dominant relative to the component p ni and thus arbitrarily strongly damped. Eventually this 
damping must reduce its kinetic energy so that it will then scale slowly (since £ — > 1) and begin to catch up again 
with the p n component. It is not surprising then to find that there are in fact, for p n ^ 0, a quite different set of 
solutions to (|[)-(0) for the exponential potential [|l{||o|, in which the energy in the scalar field mimics that of the 
dominant component, contributing a fixed fraction of the energy density determined by A given by 

fi, = -^- = £ ^ocl/a" e = l-g ( 9 ) 

P<j> + Pn A 2 6 

This solution is also an attractor pi] , with <j> again evolving logarithmically in time as given by (^) (with only 4> 
differing) . Fig [l] illustrates the evolution towards the attractor starting from initial conditions with the scalar field 
energy very dominant (obtained by a numerical evolution of the homogeneous equations of motion for the exponential 
scalar field with A = 4 and components of radiation and matter which are equal at a = 1). The scalar field energy 
first scales rapidly, approximately as 1 /a 6 as indicated by the solution (||) , until it has undershot the energy density 
in the radiation. It then turns around and starts scaling much slower than radiation or matter until it again catches 
up with them in the matter dominated regime, and then settles down at the fraction given by @ with n — 3. As 
anticipated above the main feature - the turn around in the scaling - can be understood to arise from the increase in 
damping as the radiation/matter becomes dominant. One can easily see quantitatively how this comes about: The 
energy scaling as 1/a 6 scaling occurs while the V'(<fi) term in (||) is very subdominant and we have then cj) oc 1/a 3 . If 
the dominant component on the RHS of (||) is radiation/matter one then obtains the evolution of the field 

<t>{t) = o + o io(l-(y) 2 ) Pn«P4> (n = 3) 
0(t) = o + 2<M o (l-(^)5) p n «p^ (n = 4) (10) 
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FIG. 2. The evolution of the fractional energy density in the scalar field for a selection of As 



In contrast to the case of the logarithmic dependence on time of the attractors, the potential term evolves in either 
case more slowly than the first two terms in (g) . This results in a slow-down in the scaling of the scalar energy which 
drives it back eventually towards the dominant component. We will discuss in greater detail below what determines 
the duration of this transient period in which the scalar energy is very sub-dominant. We will also see below that the 
re-entry in this example in the matter-dominated epoch is just a result of the inital conditions. For A > 2 the attractor 
exists both in the radiation and matter dominated epochs. Given that the relative contribution of the scalar field 
energy and the dominant component differs so little in the two epochs (by a factor of 3/4) one would anticipate that 
the scalar field, if established in the attractor in the radiation dominated epoch will match tightly onto its asymptotic 
value in the matter dominated epoch. That this is indeed the case can be seen from Fig |[ which shows the evolution 
of the fractional contribution to the energy density for a range of values of A (with the scalar field prepared in the 
attractor at the beginning of the evolution). 

The conclusion from this discussion is that, in the presence of matter or radiation, a scalar potential which is 
flatter than a simple exponential with A = v3 will support modes with energy scaling slower than matter, and always 
asymptotically dominate. Simple exponentials with A < 2 have attractors in the radiation and matter dominated 
epoch given by (||), and those with \/3 < A < 2 in the matter dominated epoch only. A scalar field rolling down 
a potential steeper than these simple exponential e.g. ~ l Mp will clearly always decay away asymptotically 

relative to radiation/matter. We also saw that oscillating modes scale in a way which is independent of what other 
component is present, and there are no attractor solutions. 

The attractor solutions for the exponential scalar field are therefore a special case. Before embarking on the full 
treatment of the cosmology in which the attractor (|^) describes the unperturbed limit, we discuss further why this 
attractor solution is particularly interesting from a cosmological point of view: By virtue of its being an attractor, 
it avoids the tuning problems which are inevitably a feature of all other scalar field cosmologies which have been 
proposed. First, however, we discuss briefly another attractive aspect of the required potentials. 



Scalar fields with simple exponential potentials occur in fact quite generically in certain kinds of particle physics 
theories. Because of the existence of the power-law inflationary solutions (||) such theories have been studied in some 
detail by various authors, with attention focussed on models and parameters which lead to inflation. We will not 
carry out any detailed analysis here of particular theories which would produce the exact parameters required for 
the present case, but limit ourselves to a brief review of such theories and the observation that the parameter values 
which we require are quite reasonable in the context of such theories. 



B. Theoretical Origin of Exponential Potentials 



G 



The first are Kaluza-Klein theories in which the fundamental theory has extra dimensions, which are compactificd 
to produce the four dimensional world we observe. In the effective four dimensional theory scalar fields arise which 
correspond to the dynamical degree of freedom L associated with the 'stretching' of the internal dimensions. The 
potential is generically exponential because the kinetic term has the form (L/L) 2 , while the potential 'stretching 
energy' is polynomial in L. A simple example p^JTII is a six-dimensional Einstein-Maxwell theory where the two 
extra dimensions are compactified to S 2 with radius L. Defining a field <& 

<5> = 2M P hiL (11) 

results simply in four-dimensional Einstein gravity minimally coupled to this field which then has the potential 

^($)cxexp(- ] |-)[l-exp(- ] |-)] 2 (12) 

In the regime of large $ we have V(&) ~ exp(— i.e. effectively a simple exponential potential with A = 2. 
Another example studied in detail in |5l]-|33| is the case of a pure gravitational theory with higher derivative terms in 
an arbitrary number of dimensions, D + 4, where spontaneous compactification arises because of the higher derivative 
terms. The effective four dimensional theory is in this case considerably more complicated, but again includes a scalar 
field <& defined as the logarithm of the "radius" of the compactified D-dimensions which gives, in certain regimes, a 
potential driving the dynamics of the zero mode of the field of the form of (§) with - exp(-£> ). 

Another set of models in which such potentials appear is in supergravity and superstring theories [B5| . One of the 
most studied is the Salam-Sezgin model with, N — 2 supergravity coupled to matter in six dimensions. It predicts 
PH the existence of two scalar fields $ and T with a potential of the form 

V(<S>,?) = V(T)eM-^^) (13) 

For T 3> it corresponds to a potential of the required form we want, with A = \/2. Further examples relevant to 
inflation are given by the authors of fl36[ , who find two scalar fields with A = V2 and A = in and N — 2 model, as 
well as in an N = 1 ten-dimensional model with gaugino condensation ]j37| . 

A further class of theories in which such fields arise is in higher order gravity. There is a conformal equivalence (see 
p8| , and references therein) between pure gravity described by a lagrangian which is an analytic function of the scalar 
curvature R and general relativity plus a scalar field with a determined potential, and in the presence of matter the 
metric associated with the latter description is the physical one 63]. For example f3§| l in d dimensions the potential 
in the case of a lagrangian quadratic in R is 

m ocex P (^i^)[l-cxp(-^^)P (14) 
and, for d — 4 and a polynomial lagrangian of the form J^ w _ 1 a n R n , it is 

V@) = A, cx P (-2^) + A n exp(-2 J-)[l - exp(^)]" (15) 

where the A n are rational coefficients determined by the a n . As $ — » +oo we have V ~ exp[(fc — 2)<E>/Mp]. 

All previous analysis of the cosmological consequences of the existence of such scalar fields has focussed on inflation, 
which is realized in the case of the simple exponential for A < y/2. The fact that many of these models can at best 
give A < y2 rather than the considerably flatter potentials needed for inflation (to naturally give a large number 
of e-foldings and a nearly flat spectrum of perturbations) meant that these theories provided a general motivation 
rather than a realistic model. In |3(|, for example, an extra ad hoc damping term was introduced to produce a more 
satisfactory model from supergravity motivated models. In the present context this is not the case: We will see that the 
most interesting range for structure formation is € [0.08, 0.12] in the matter era which approximately corresponds 
to the range A € [5,6]. Although these values are a little larger than in the simplest models we have reviewed, 
some of the models above clearly can lead to parameters in this range, and, for example, the number of compactified 
dimensions required to give these values in the first type of theory we reviewed is certainly not unreasonable in the 
context of superstring theory. In terms of the theoretical origin of the potential we study therefore, it is not just the 
form of the potential which is one which arises naturally, but it may also be that the precise value of the single free 
parameter in this potential is a natural one. Clearly further analysis of such models would be required to make a 
stronger statement than this, and to see if the cosmological effect of these fields which we study in this paper might 
ultimately be used to give us specific hints about fundamental particle physics. 
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C. Nucleosynthesis Constraints 



In considering the possibility that some significant fraction of the energy density of the universe may be in a 
homogeneous mode of a scalar field, the earliest constraint comes from nucleosynthesis. Such a contribution from a 
weakly coupled scalar field enters in determining the outcome of nucleosynthesis (the primordial densities of the light 
nuclei) only through the expansion rate. Adding a component increases the expansion rate at a given temperature. 
The dominant effect of such a change is in its effect on the ratio of neutrons to protons when the weak interactions 
keeping them in equilibrium freeze-out at a temperature of ~ lMeV. The range of expansion rates at this temperature 
compatible with observations is usually translated into a bound on the number of effective relativistic degrees of 
freedom at ~ lMeV. Taking AN e ff to be the maximum number of such degrees of freedom additional to those of 
the standard model (with three massless neutrinos) the equivalent bound on the contribution from a scalar field is 

(where 10.75 is number of effective degrees of freedom in the standard model). A wide range of values for the upper 
bounds on AN e f / exists in the literature, the variation being due both to the data taken into account and the methods 
of analysis used. We will not attempt here to review all the issues involved and accept one or other bound as definitive 
(see , for example, 0] for a full discussion). The tendency in the last few years has been towards less restrictive 
bounds than were generally thought correct previously. A typical value now used is a bound of AN eft — 0.9 which 
is given by various authors |^0|, or even a more conservative one of AN e ff — 1.5 by others pd| , |42]| . This range 
corresponds here to 

Q^lMeV) < 0.13 - 0.2 (17) 

This is the range of values we will take when we discuss nucleosynthesis in the rest of the paper. The result for more 
restrictive nucleosynthesis analyses can be read off from ([ji]) . Note also that we have assumed here that the standard 
model number of relativistic degrees of freedom at lMeV. The strict experimental lower bound on this number is in 
fact 9, given that the upper bound on the mass of the r neutrino is 18 MeV ]2q] . In the case that the r neutrino is 
non-relativistic (and decays before nucleosynthesis) the bound ([l7]) is changed to Q^ilMeV) < 0.27— 0.33. 



D. Scalar Fields at Late Time and Fine- Tuning 

In modifications of standard inflationary flat cosmologies |l6|-|lS|] involving a contribution from a scalar field which 
have been considered to date, attention has been focussed exclusively on the case that the scalar field contributes 
significantly only at recent epochs, at the earliest well after the transition to matter domination. The main reason 
for this is that one of the motivations for many models has been to produce a contribution at late times which scales 
slower than matter and dominates asymptotically, producing effects very similar to that of a cosmological constant. In 
this case there is unavoidably the same sort of tuning as involved in the cosmological constant model: One requires an 
initial energy density in the scalar field which is characterised by the energy density in the universe at recent epochs. 

Another kind of model which has been considered is that in which the energy density in the scalar field scales like 
matter asymptotically, implemented in an oscillating mode of a sinusoidal potential. If the field lies initially away 
from the minimum, it becomes important once the curvature of the potential is comparable to the expansion rate and 
initially behave like a cosmological constant before rolling down the potential and scaling asymptotically like matter. 
This model also involves the same sort of tuning, since a small energy scale must be introduced to single out this 
late time at which the scalar field becomes dynamically important. It has been argued in |L2] that there are particle 
physics motivations for the introduction of such a potential with such a characteristic scale, and that there is nothing 
'unnatural' about models of this type. 

The special case we have discussed of the attractor in the exponential potential is quite different, simply because it 
is an attractor. In contrast the only attractor in the case of a field asymptotically scaling slower than matter is the 
pure scalar field dominated cosmology, and the tuning referred to is just that required to make the coupled system 
of scalar field and matter approach this attractor just at the present epoch. For the models with an oscillating mode 
there is no attractor, and the tuning consists in fixing the energy density in the scalar field to be comparable to that 
in matter when it starts scaling like matter. With the attractor for matter and radiation in the present case, the 
scalar field plus matter/radiation will always end up in this solution asymptotically. In this paper we will take this 
solution to apply from the beginning of our simulation of structure formation, deep in the radiation era at a red-shift 
of z ~ 10 7 when the temperature is ~ 100eV^. We now address the question as to what range of initial conditions 
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on the scalar field in the early universe will give rise to this behaviour. In the course of this analysis we will also 
determine what sorts of initial conditions are compatible with the late entry (in the matter era at a red-shift of ~ 70) 
to the attractor discussed in 



E. Initial Conditions in the Early Universe and the 'Self- Tuning' Scalar Field 



The assessment of what are 'natural' initial conditions for a scalar field in the early universe requires of course a 
particular framework within which to address the question. What we need to determine is essentially just the kinetic 
and potential energy in the scalar field at some early time relative to that in radiation (and therefore matter) . Given 
that we are working in the context of inflation-motivated flat homogeneous cosmologies (and will assume Gaussian 
adiabatic perturbations of the type produced by inflation) we assess the question within this framework. The energy 
density in radiation is then determined by how the universe is reheated after inflation. We consider both the usual 
scenario for reheating by decay of the inflaton and then an alternative scenario introduced in J47| . The reason for our 
detailed analysis of this second non-standard case will become clear below. 

First consider the standard reheating scenario. We suppose there is an inflaton field and the scalar field with 
an exponential potential V e^ x ^/ Mp with A > 2 (i.e. with the attractor in the radiation/matter epochs which we 
will consider). Let us consider a typical inflationary model e.g. chaotic inflation. The simplest and most natural 
assumption for the relative energy densities at the onset of inflation is 

\tf ~ V(<f>) < V mf (18) 

where Vi n f is the energy density in the inflaton. We make the latter assumption since it is required for the onset of 
inflation. The dynamics of the two fields are then described by the equations (g)-(Q), but with p n now the total energy 
of the inflaton and n a function of the inflaton n = 6(1 — £in/)> where is the ratio of the potential energy of the 
inflaton to its total energy. For the same reasons as in the case discussed above when p n describes matter or radiation, 
there is an attractor solution given by (^|) with n = 6(1 — (assuming n changes slowly which is the case for 

slow- roll inflation) . The energy in the scalar field scales so slowly because its roll is strongly damped by the inflaton, 
with (||) specifying the exact ratio of energies at which the damping slows the roll to give precisely the same scaling 
as for the inflaton. In inflation driven by an exactly constant energy density, this ratio is zero and just represents the 
asymptotically approached solution in which the scalar field rolls away to an arbitrarily large values of the field after 
an arbitrarily long time. In any realistic model of inflation however the inflaton must roll in a non-trivial potential 
in order to exit from inflation and £i n f ^ 1. For example, in chaotic inflation in a potential ~ <f>^ n t, one has, in the 
slow-roll inflationary regime, n = |(^^) 2 - Once inflation commences, say at 4> in f — 4>1 n fi the ener gy m the scalar 
field will be driven towards 

"♦"A'" 9 V? n/ j ~9iV e (19) 

where fi™ is the fraction of the energy density in the attractor in the matter era, and iV e w i s the number of 

e-foldings of inflation. As the inflaton rolls down the potential n increases and the fraction of energy in the attractor 
grows. 

Starting from initial conditions like ( p"8| ) the roll of the scalar field is very rapidly over-damped due to the rapid 
red-shifting of its initial energy density. It will always then scale much slower than it would in the absence of the 
inflaton, with In the next section we will see that this is enough to ensure that its energy density relative to the 

inflaton will never drop substantially below that in the attractor. (The transient with sub-dominance we observe in 
that case results from the field having initially evolved without the p n component playing any role for a long period). 
Without analysing in detail how precisely the energy density can track the (growing) attractor value during inflation, 
we thus conclude that the fraction of the energy in the scalar field at the end of inflation will be bounded below by 
Sl% as given by ©. 

From Fig || (and our previous discussion of Fig |]) we see that this conclusion is enough to establish that, after 
inflation, the attractor for the scalar field in the radiation dominated epoch is approached within a few expansion 
times of the end of inflation. When the inflaton starts oscillating and/or decays it begins to scale like matter or 
radiation, while the scalar field remains overdamped and its energy almost constant until the attractor is attained. 
In Fig [3] we assume an abrupt transition from inflation with essentially constant energy to radiation type scaling at 
scale factor a = 1, and take the initial fraction of the energy in the scalar field to be 3 x 10 -4 . The energy density 
approaches that in the attractor by a ~ 10, and then oscillates about it as it approaches it asymptotically. Without 








lo ^ a > log(a) 
FIG. 3. In the left panel we plot the evolution of the energy density in the scalar field (p$) with A = 4 and in the remaining 
matter as a function of scale factor (the units are arbitrary), in the case that the scalar field is initially sub-dominant. In the 
right panel we plot the evolution of the fractional density in the scalar field. 

further detailed analysis it is clear that for a typical inflationary model, in which the reheat temperature is at most 
a few orders of magnitude below the GUT scale, the attractor will be approached to great accuracy not only by the 
time scales relevant to structure formation, as we assume below, but also long before nucleosynthesis. 

Late entry into the attractor for the scalar field, in particular late in the matter dominated epoch as assumed in 
the case of this type treated in ]l8| ] , is therefore certainly not what one would expect from typical inflationary models 
with the usual method of reheating. Inflation is however not embodied in a specific model or set of models, and it 
would certainly be possible (and easy to see how) to devise a model in which the energy density in the scalar field at 
the end of inflation could be tuned to any required value. 

Given that the attractor solution is typically established by nucleosynthesis, one must satisfy the constraint discussed 
above. Converting ([l7]) to one on the contribution from the scalar field at late times (in the matter dominated epoch) 
givesQ 

^ < 0.1 -0.15 or A > 5.5 -4.5 (20) 

It is because this quantity would seem to be too small to play any important cosmological role that this kind of 
attractor solution has been disregarded by the authors who have noted these solutions^. That this is not in fact the 
case we will see in greater detail in the later part of this paper. Because the contribution is important for a long time 
it can have just as significant an effect as a larger contribution which plays a role only at late times. 

The initial motivation for our detailed study of this model was not, however, the realization that the attractor 
might consistently be established by nucleosynthesis for cosmologically relevant cases. Rather our starting point was 
the observation that there are interesting and viable cosmologies in which scalar field energy in a rapidly scaling mode 
could dominate over radiation in the very early universe, and that in the case that this field is an exponential of the 



2 The upper bound here which corresponds to the analysis of jllj is actually given explicitly in this form (for the exponential 
attractor) in p2| . 

3 Ratra and Peebles |jl6l give this explicitly as the reason for the attractor solution model being "phenomenologically untenable" 
(p. 3416). Wetterich in |20| considers the homogeneous cosmology of the case in which the exponential field dominates at late 
times; in order to satisfy the nucleosynthesis constraint he discusses the possibility of late entry brought about by a tuning of 
the initial scalar energy density to an appropriately small value, and also the possibility that A, rather than being constant, 
decreases between nucleosynthesis and late times. 
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type relevant to late time cosmology, there may be, as shown in Fig [j] a transient period between the two epochs (of 
domination by the scalar field, and the late time attractor) lasting many expansion times in which the scalar field 
energy is negligible. If this transient period includes nucleosynthesis the constraint (|2(]) would not apply. We now 
discuss this model and examine how the time of entry into the attractor depends on the parameters in the model. 



An epoch dominated by a scalar field in a mode scaling faster than radiation (or kination |44[])) comes about by 
construction in an alternative theory of reheating suggested in |47|j . Instead of rolling down to the minimum of a 
potential, oscillating and decaying, as envisaged in the standard reheating scenario, the inflaton field can roll into a 
steeper potential supporting a mode scaling faster than radiation (i.e. an exponential with A > 2 or steeper). Rather 
than being, as is often stated, completely 'cold and empty' at the end of inflation, the universe contains a component 
of radiation created simply by the expansion of the universe (with energy density ~ H 4 ). Although initially very 
sub-dominant relative to the energy p<f, in the scalar field (H 4 / p^ ~ p^/Mp) a transition to a radiation dominated 
cosmology will take place at some subsequent time since the radiation red-shifts away slower than the energy in the 
scalar field. 

We restrict ourselves to the case that the relevant field is the simple exponential with A > 2 (i.e. £1^ < 0.75). We 
evolve the system of scalar field plus radiation and matter forward in time from the end of inflation, at which time 
we take the initial conditions to be specified by 



where Hi is the expansion rate at the end of inflation. The initial condition assumes an abrupt end to slow-roll 
inflation i.e. when the inflaton begins rolling in the region in which the potential is exponential, the potential energy 
is still dominant. The choice e = 1CP 3 corresponds to the simplest estimate of the initial radiation density, taking 
it to be dominated by the radiation at the horizon scale at the end of inflation (with 'temperature' T ~ H/2ir fl46| , 
see ^,^5[ for a more detailed discussion). We will see below that the results we are interested in here are not very 
sensitive to these choices. 

In this model with the exponential potential there are then just two parameters: Hi, the expansion rate at the end 
of inflation, and A (or, equivalently, fi^ = 3/ A 2 in the attractor in matter domination). These entirely specify the 
post-inflationary homogeneous cosmology, which evolves, as illustrated in Fig [l], from the scalar field dominated phase 
through a transient into the late time attractor. What we want to determine here is the time at which the attractor is 
approached as a function of Hi, and the range of parameters for which the model is compatible with nucleosynthesis. 

The results, which are summarized in Fig |^, can be understood after some closer examination of the solutions. The 
first feature to note is that the nucleosynthesis constraint, denoted by the hatched region, allows a range of Hi at 
given fi^ which is (i) only bounded below for < 0.15, (ii) bounded above and below for fta > 0.15 in a range 
which pinches off as f2^ — > 0.5, so that for larger values of there is no parameter space in the model compatible 
with nucleosynthesis. 

This can be understood as follows. The lower bound comes from the requirement that the initial regime of scalar 
field dominance end sufficiently long before nucleosynthesis. In its initial phase the scalar field approaches rapidly a 
mode in which it scales approximately as 1/a 6 for < 0.5 (or as given by (^) for f2^ < 0.5 (A > V&))- In order that 
the radiation p l rad come to dominate by the time its temperature (T ~ Hiai/a) corresponds to that at nucleosynthesis 
(T ns ~ 1 MeV), this means (from @) that 



where T enc i is defined as the temperature at which the scalar field energy becomes equal to that in the radiation, 
and we took T en d = 3 MeV to obtain the numerical value. The higher value of the lower bound in Fig |] results 
from the fact, since our initial conditions arc still inflationary, the scalar field takes a short time to attain the rapidly 
scaling behaviour assumed in deriving (p2|). The increase of the lower bound with is explained in the same way - 
the transition to the 1/a 6 scaling takes longer as the potential flattens (as A decreases). From (^2|) we also can see 
quantitatively the weak dependence on e. 

As we observed in Fig |l| the scalar field dominance is followed by a transient period in which it is very sub-dominant, 
before approaching the late-time attractor. The upper bound on Hi which appears at £1$ = 0.15 is just the requirement 
that, for > 0.15 (using the weaker limit of the bound in (p"7j)), this period lasts until after nucleosynthesis. This 
gives an upper bound because the transient ends earlier according as the first phase of scalar field dominance does. In 



F. Late Entry in an Alternative Model of Reheating 



V(4, )=3M P H? O = O p I3A = eH? e = HT 3 



(21) 




(22) 
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0.2 0.4 0.6 



FIG. 4. Reheating by kination in a simple exponential potential: The solid region is that excluded by nucleosynthesis 
constraints. The solid line (dash-dot, dotted) show the models for which the attractor is established at the beginning of 
structure formation (matter domination, today). 
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FIG. 5. The fractional energy density as a function of temperature for Q<p = 0.3 and three values of Hi. 




Fig we see the evolution of the fraction of the energy density as a function of temperature for various different Hi . 
For Hi = 10 12 GeV the nucleosynthesis bound is violated because we enter the attractor too early ; for Hi = lO s GeV 
it is violated because we exit the first phase too late. Why the allowed range of Hi 'pinches off' as we go to = 0.5 
is that the duration of the transient decreases as we approach this value. This can be seen clearly from Fig [6] which 
show, for £1 = 0.1 (left panel) and £1 = 0.5 (right panel) the temperature T en d (when the scalar field energy first 
equals that in the radiation) and at T att , when the attractor is established. This behaviour can be understood by 
looking back to the attractor solutions (||). In the late time attractor the parameter £ = V(4>)/ p<j, must reach a certain 
determined value; as indicated by (JsJ) the value A = \f& separates two quite distinct regions in which £ behaves quite 
differently. For A < y/6 the p n = attractor tends to a finite value; in the steeper potential with A > s/& the potential 
energy 'runs away' faster and £ decreases rapidly towards zero. As long as £ ~ the scalar energy continues to scale 
away as 1/a 6 , and it may take a long time for the kinetic energy to catch up again on the potential energy (which 
essentially stops decreasing once the field's roll is damped by the radiation or matter). The large 'undershoot' and 
long period of sub-dominance of the scalar energy thus results from the fact that the exponential field with A > \/6 
evolves for a long time damped only by its own energy density, which allows the field to run away so fast that the 
potential energy decreases enormously. On the other hand for £1 = 0.6 (A = \/E) the scalar field in the first phase 
follows the p n = attractor (@) with £ = 5/6, and only a little adjustment is needed to find the late time radiation 
dominated attractor with £ = 2/3. The nucleosynthesis constraint is therefore never satisfied in this case. 

Also shown in Fig ^ are the curves corresponding to 'entry' to the attractor at different temperatures - at lOOeV 
(solid line), at matter-radiation equality (dash-dot line) and today (dotted line). The undashed region lying to the 
right of the first of these lines corresponds to the parameter space (for this model) we will consider in the rest of this 
paper, which is consistent with nucleosynthesis and with our assumption of the validity of the attractor with (^|) as 
our homogeneous solution. The line corresponding to re-entry at nucleosynthesis is not given, but one can infer that 
it would mark out a substantial part of this region (bounded by the solid line), which corresponds to entry into the 
the attractor after nucleosynthesis but prior to T = lOOe^. There is also a region, bounded by the first and third 
line, which we do not describe and which might be of relevance to late time cosmology. The region to the left of the 
dotted line corresponds to regions where the scalar field would still be irrelevant today, and therefore the existence of 
the exponential field has no implications for late time cosmology. 

One can also determine from Fig |]that, for entry in the matter dominated epoch, the largest possible value of Q.^ 
in this model is approximately 0.3. This places an upper bound on what contribution to the energy density might be 
expected to come from this type of solution, if one is not to use the sort of tuning which, as we have emphasised, the 
model has the merit of being able to avoid. (In particular we note that we cannot obtain the case = 0.6 of |lq| .) 

Is there a 'natural' value for Hi within the allowed range? Hi is related to the inflaton energy at the end of inflation 
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by Pl ~ (HM P ) 2 , and the range H t G [10 8 , 10 16 ]GeV corresponds approximately to p] /4 £ [10 13 , 10 17 ]GeV. With a 
more developed model with a full ansatz for the part of the potential which supports inflation, one can relate H t to 
the amplitude of density perturbations, and require that they match those observed in the CMBR by COBE. The 
simplest calculation (for a very flat inflaton potential) would give p]^ 4 ~ 10 16 GeV, or Hi ~ 10 14 GeV. However the 
result is not model independent. In [ ft5| a particular model is constructed in which H can be as low as the value at 
which the phase of scalar dominance ends just before nucleosynthesis. This parameter range for the model has the 
particular interest that it leads to quite radical modifications of physics at the electroweak scale, since the expansion 
rate can be up to five order of magnitude greater than in the radiation dominated case at T ~ lOOGeV '. 

In the rest of this paper we study only cosmology starting from T ~ lOOeV^. We will assume the existence of the 
simple exponential, and take the homogenous cosmology about which we perturb to be given exactly by the attractor 
at this initial time. A is thus the sole adjustable parameter additional to those of SCDM. To allow the simplest 
direct comparison of structure formation in our model with SCDM we take the fiducial value VL h 2 = 0.0125 for the 
baryon fraction. In the present study we will not consider the effect of varying this parameter. If we were to do so 
we would have to be more definite about our assumptions about nucleosynthesis, since the allowed range of f2f,/i 2 (or 
equivalently baryon to photon ratio rj) depends on the expansion rate at nucleosynthesis. The canonically quoted 
range, which corresponds to that consistent with observations for the case N v = 3, narrows as the expansion rate at 
freeze-out increases (since the fraction of Helium produced increases as the expansion rate does). As the expansion 
rate reaches its upper bound, the rj required for consistency with observations is pushed towards its lower bound. 
Having used the constraints of structure formation to determine the best value of fi^ in our model, we will comment 
again in our conclusions on the consistency of entry to the attractor prior to nucleosynthesis. 



III. EVOLUTION OF PERTURBATIONS 



We now analyse the evolution of perturbations in this cosmology. The structure of this section is as follows. We 
firstly write down the equations for the unperturbed cosmology in conformal coordinates which we use in the rest of 
the paper. For complet eness , in III B we present the equations for the evolution of cosmological perturbations in the 
synchronous gauge. In III C we derive analytic solutions to the evolution of perturbations on sub- and superhorizon 
scales in the radiation and matter era. It becomes apparent that they are very similar to those in a cosmology in 
which part of the matter content is in a massive neutrino (M DM). It i s thus instructive to identify the key differences 
between our scenario and the MDM model, and we do so in [ill Dj . In |III Fj we look at the evolution of perturbations 
in the radiation and how it affects the temperature anisotropy angular power spectrum. 



A. Equations of Motion and Initial Conditions for the Background 

Normalising the scalar field in units of Mp and defining the origin of <& so that the scalar potential is Mj,e~ A *, 
the equation of motion for the unperturbed background scalar field is, in conformal coordinates for which ds 2 = 
a 2 [—dr 2 + 8i J dx t dx ] ], 

4> + 2Hj) - a 2 XMp exp(-A^) = 

V 2 = ^a 2 (p^ + PM ) (23) 

where a is the scale factor, Ti is the conformal Hubble factor, dots denote derivatives with respect to conformal time 
t, pm is the energy density in matter and radiation, and 

1 r 2 +M 2 exp(-A0) (24) 



P<p 



Mi 



2a 2 



In the attractor (ffl we have 



PK 



T2 Ptotal = 3Mp I — I (25) 

p<t> 6 ptotal a \ a J 

where px — -^z4> 2 is the kinetic energy density of cf>. We use these to fix the initial conditions on <fi and </> to be 



H 



= — x lo g 



/4n(6 - n)H 2 



a 2 X 2 



(26) 



where n = 4, since we begin the evolution deep in the radiation dominated epoch. 
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B. Linear Perturbation Theory 

We use the notation and results of Ma & Bertschinger 0] , with the modifications brought about by the addition 
of a scalar field; we present the full set of equations but we refer the reader to [B and |^] for details on how to solve 
them. The formalism we work in is the synchronous gauge where ds 2 = a 2 [—dr 2 + (5ij + hij)dx l dx : >]. Restricting 
ourselves to scalar perturbations, the metric perturbations can be parametrised as 

ha = J dPke^fi&hQL, t) + (k 2 k, - ify )6»j(k, r)]. (27) 

Pressureless matter (cold dark matter) has only one non-zero component of the perturbed energy-momentum tensor: 

£T ° = - Pc 5 c (28) 

and it evolves as 

So = -\h (29) 

Radiation can be characterised in terms of a temperature brightness function, Ay(k, n, t) and polarization brightness 
Ap(k, fi, r). These brightness functions can be expanded in Legendre polynomials of k • n: 

oo 

A T (k, fi, r) = Ysi-iYW + 1)At/(A, r)Pi(k ■ fi) 

l=o 
oo 

A P (k, fi, r) = ^H) £ (2£ + l)A Pe (k, r)P e (k ■ fi). (30) 



Defining the density, velocity and shear perturbations by 

Sry = ApO Of = 

the perturbed energy-momentum tensor for radiation is 



<5 7 = A T0 <9 7 = -kA Tl u 1 = ^A T2 (31) 



STq = —PfS. 



7 



.0 4 



ik l 5T? = -p 7 7 
ST* = \ Pl 6 7 + Ej 



(k,-kj- - -Sij)E) = --p 7 cr 7 (32) 



3 

1. 4 
3 

Thomson scattering couples the radiation and baryons, and the latter have a perturbed energy momentum tensor: 

ST" = -p b 5 b 



ik l 5T® = ^ Pl 9 b (33) 



The evolution equations for radiation are: 

4 2 • 

(5-v = Ory — fl 

7 3 7 3 



# 7 = k 2 (-gfi-y ~ &~^J an e (TT(0b — &y) 

8 3 4 8 

2 ^ = 15^-5 fcA73+ 15 /l+ 5'' 

9 1 
--an e er T er 7 + -an e ir T (Apo - A P2 ) 

k 

A T £ = 2 £ + i ^ A T(e-i) - (t+ l)A T (i+i)] - an e a T A T e 
^r^J [^P(/-l] " (* + 1)Ap(£+i)] 

+fln e ir[Apf + i(A T2 + A P0 + A P2 )(<5<m + %)] (34) 
2 5 



A 
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and for baryons 



h = -6b - 7^h 

6 b = -W b + Ran e a T {6 1 -e b ) (35) 

ut is the Thomson scattering cross section, n e is the electron density, and R — The evolution equations for 

massless neutrinos can be obtained from those for radiation by setting Ap = R = a? = and T — > v. 
The perturbation ip in the scalar field about the homogeneous solution <f> has the equation of motion 

<p + m<p-'V 2 <p + a 2 V"<p + -<j)h = (36) 

and gives rise to perturbations in the energy momentum tensor which are 

a 2 5T° = ~4>ip ~ a 2 V'<p 
-a 2 di5T° = <j)V 2 (p 

a 2 8Tl = 3<^> - 3o 2 V"V- ( 37 ) 

Finally we have the perturbed Einstein equations 

k 2 V - ^Hh = 4nGa 2 6T° 

k 2 fj = 47rGa 2 zki^° 
h + 2Hh - 2k 2 -q = -8nGa 2 5Ti 

h + 6ij + 2H(h + 6fj) - 2k 2 rj = +24vrGa 2 (k t k J - -8 i3 r)Sj (38) 



Having defined the evolution equations for the perturbations, all that remains to be specified are the initial condi- 
tions. We are working within the context of inflation and consider an initial set of perturbations of the type generically 
predicted by the simplest such scenarios: We assume a Gaussian set of adiabatic perturbations with a scale invariant 
power spectrum. This completely defines the statistical properties of the ensemble. Putting the adiabatic perturba- 
tions in the superhorizon growing mode gives [Q : 

5 1 = -^C(kr) 2 

S c = --h = S b = |<?„ = ^(5 7 

7 23 + 4i?„ 18 V ; 

- - -(far 



3(15 + 4R„) 

V = 2C - J + *f C(krf (39) 
6(15 + 4R V ) 



and, as we will see in more detail below, 



5A 

5 C = —ip (40) 
(and C is the overall normalisation). All the remaining perturbation variables are set to zero initially. 

C. Asymptotic solutions to the evolution of density perturbations 

We can arrive at a simple understanding of how perturbations evolve in this scenario by considering the simplified 
case in which ther is only cold dark matter, radiation and the scalar field. The evolution equations are then: 
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6 C + H6 C - -H 2 {n c Sc + 2Q r 5 r ) -2<j)tp + a 2 V'ip = (41) 

S r + \k 2 5 r - p c = (42) 
ip + 2H(p + k 2 ip + a 2 V"tp - #c = (43) 

where fix = px/ptot (X = C,7, f, <j>). 

Consider first the superhorizon evolution. In this limit we set k 2 = and (assuming adiabatic initial conditions as 
above) 5 C = |<5 7 . Using the scaling solutions for the homogeneous mode of the scalar field, in the radiation era we 
have: 

o c + -o c o o c - —ip - -T—ztp = 

T T z AT AT 

2 4 4 • 

cS + -y>+ -99- -S c = = (44) 

T T A 



9 _n -lzbV-15 + 64/A^ 

These equations are solvable and give o c and y as a linear combination of r , r and r 2 . Taking only 

the coefficient of the growing mode to be non-zero, we find (as stated above) S c — ^ip = At 2 . On superhorizon scales 
in the matter era we have 

•• 2- 6(1-4), 12 18 

S c + -o c 5 Sc - t— <i> - "TT^ = 

t t z At at z 

•• 4 • 18 6 ■ 

tp + -ip + —ip - —6 C = = (45) 

T T AT 

9 _Q -6± \/-7+12/A 2 

which can be solved to give o c and ip as a linear combination of r , r and r 4 . With the assumption 
therefore of initial adiabatic perturbations in the growing mode in the radiation era, we have 5 C = 2 ^ L <-P oc t 2 . We 
conclude, therefore, that the evolution of perturbations on superhorizon scales are exactly the same as in a scalar 
field free universe. This is just a manifestation of the fact that, on very large scales, the dominant interaction is 
gravitational, which is blind to matter type. 

Next we turn to the subhorizon evolution of perturbations, i.e. fer> 1. Consider first the radiation era. Assuming 
that the gravitational feedback is unimportant on small scales, we have 



1. s 4 

1 



5 C + -S c - — — ^- L S 1 - -^-(p - -^-<p = 

T T z AT AT 



S 7 + ^fc 2 <5 7 ~ 
2 

(p + -ip + k 2 <p ~ = (46) 

T 

The last two equations are easy to solve, giving <5 7 oc e ± ^ T and ip oc -^Ji(fcr), -^Ni(kT) where J M and are 
spherical Bessel functions. Clearly these have little effect on 5 C and one can drop these terms from the first equation 
to get 8 C oc Ci, logr. The subhorizon evolution in the radiation era is thus essentially the same in a universe with 
and without a scalar field. 

In the matter era on sub-horizon scales the equations reduce to 



T T z AT AT 

6 



c23 H — cp + k 2 ip - --S c = (48) 

T AT 



To a first approximation we can discard the gravitational feedback term in Eq. 48 . The solutions to this wave equation 



are ip oc ^372 J| (fer) , -^j^N^ikT) i.e. oscillatory solutions with decaying amplitudes. Clearly this will contribute little 
to the growing mode and we can drop the last two terms in Eq. [l7]. The equation is easy to solve and, using = 

we find the growing mode solution 5 C oc r 2+ 2 25 s1 *). The subhorizon growing mode is therefore suppressed 

relative to that in a scalar field free universe. 

Having derived these approximate forms we now have a rough idea of what the mass variance per unit in In k is, 
where this is defined as 
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FIG. 7. On the left panel we plot the evolution of S c for two wavenumbers in a sCDM cosmology (dashed line) and in a 
universe with = 0.1. In the right panel we plot n e ff for k — O.lMpc" 1 for the same cosmology. The upper (lower) dotted 
line is the asymptotic matter era solution in the sCDM (S7^ = 0.1) case. 

A 2 (fc)~(|<5 c (fc)| 2 ). (49) 



One has (denning 2e = 5(-l + x /l - 

fc 4 if fc < 2zl 

i — T ° 

A 2 (fc)oc { fc 4 " 2 ' iik <rf q (50) 

constant + In fc if fc > — 

It is instructive to see the evolution of density perturbations for a few wavenumbers, one that comes into the horizon 
at around radiation-matter equality (fc = O.lMpc ) and one that comes in during the radiation era (fc = l.OMpc -1 ) 
for = 0.1. They are compared in Fig ^ to density perturbations in a universe with no scalar field and it is clear 
that there is a suppression of growth after horizon crossing. Another useful quantity to plot is the dimensionless 
growth rate n e ff — j-5 c . For a scalar field free universe we know that deep in the radiation or matter era n e ff = 2. 
For the scalar field cosmology we have deep in the radiation era n e ff — 2 but deep in the matter era n e tf = 2 — e. In 
Fig we plot n e ff for fc = O.lMpc -1 and can clearly see the different asymptotic regimes. Note that there is a long 
transient to the matter era solution in both the cosmologies considered. 



D. A comparison with the evolution of perturbations with mixed dark matter 

By analysing the simplified system we have been able to get an idea of what to expect from solving the full set of 
perturbation equations. We found that the presence of 4> in the matter era suppressed the growth of perturbations by 
an easily determined factor. The behaviour of this system is very much like that of an MDM cosmology. There one has, 
in addition to pressureless matter, radiation, baryons and massless neutrinos, a component of matter in two species 
of massive neutrinos with masses, m v of order a few eV. The evolution of the energy density is similar to that of the 
scalar field in the present case. Deep in the radiation era, these neutrinos behave as a massless species and therefore 
the energy density scales as radiation (cx 1/a 4 ), while deep in the matter era they behave as non- relativistic matter 
which scales as 1/a 3 . Unlike the scalar field case, where the transition between the two regimes is determined by the 
transition from radiation to matter domination the transition in the case of the neutrinos occurs when 3fcBT„ ~ m v , 
where fc^ is the Boltzmann constant and Tv is the massive neutrino temperature. This corresponds to redshift 
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FIG. 8. The evolution of the fractional energy density in <j> (dashed line) for a h = 0.5 universe with — 0.1, Q c = 0.85 and 
Qb = 0.05 and the fractional energy denisty in two species of massive neutrinos (solid line) with f2„ = 0.1 (massive), £l c = 0.85 



z ~ 1.8 x 10 5 . In Fig || we plot the evolution of the energy densities in the extra component in each of the two 
cosmologies. The energy density in the scalar field <f> follows the transition very tightly while the energy density in 
neutrinos becomes non-relativistic after radiation-matter equality. In the latter case there is a short period of time 
after equality when the massive neutrinos contribute less to the energy density than they do asymptotically, and as 
a result the pressureless matter clumps more strongly for this period. 

The correction to the exponent of the growing mode in the matter era is also very similar in the two models. 

In |Q it was shown that the correction for MDM is e„ = |(— 1 + yjl — ||f2„), exactly the same as we have just 
derived for the scalar field cosmology. The reason is just that the sole assumption in the derivation of this result 
was that the exotic form of matter (the <fi field in our case and the massive neutrinos in the MDM case) does not 
cluster below a certain scale. There is however again a small but important difference. For <p, the scale below which 
a given mode does not cluster is the horizon, i.e. cx ^. Once it comes into the horizon it never clusters. For the 
massive neutrino perturbations this is not the case. The free-streaming scale, when the neutrinos are non-relativistic 
is kf s = 8a 1 / 2 (m l/ /10el/)/iMpc _1 , which grows with time. Perturbations of a given wave number k are damped while 
k > kf s , but as soon as kf s becomes small enough they behave like perturbations in pressureless matter and grow. 
So perturbations in an intermediate range of wavenumbers, despite suppression during a finite period of time, catch 
up on the pressureless matter perturbations and contribute to gravitational collapse. This effect can be seen when 
we compare n e ff for both these scenarios for k — IMpc (i.e. a mode that comes into the horizon at radiation-matter 
equality). In Fig ^ we see that, while in the 0CDM case n e ff — > e, in the MDM case, perturbations in massive 
neutrinos start to grow as pressureless matter and n e ff — > 2. 

A further subtle difference between the two scenarios can be seen at the transition between super- and subhorizon 
behaviour in the radiation era. As shown above, the subhorizon evolution of S c is cx constant, lnr. The initial 
amplitude of these solutions is set by the amplitude at horizon crossing, i.e. when fcr ~ 1. In the MDM case, 
the analysis is also simple. The neutrinos are effectively massless and therefore behave exactly like radiation. So 
again the subhorizon evolution of S c is oc constant, lnr. However the transition to subhorizon evolution happens at 
smaller scales than for the tp. Indeed from the equations we see that the transition should happen when c s kr ~ 1 
(where c 2 s = i). A simple way of seeing this is by looking at the source term for the S c in the radiation era, 

AS(k,r) — S c + -5 C = d T {rn e ff)^%. On superhorizon scales S(k,r) = 1, and the time when this quantity starts to 
deviate from a constant indicates the transition from super- to subhorizon behaviour. We have 



and fi 6 = 0.05. 




(51) 
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FIG. 9. We plot the evolution of n e ff for k — O.lMpc in a MDM cosmology (dashed line) and in a universe with fi^j =0.1 
(solid line). The dotted line is the asymptotic matter era solution, e. 

in the MDM scenario where we have grouped the massless and massive neutrinos together, and we have 

s(k,r) = fain^ + nj^ + ^-ip + ^tp (52) 

4 AT It 

in the 4>CDM case. In Fig. [lO] we look at the time evolution of a mode with k = 5h Mpc -1 and see that transition 
to subhorizon evolution occurs earlier in the scalar field scenario than in the MDM scenario. As a result 6 C stops 
growing earlier and there is an additional small, but non-ncgligablc suppression factor in the radiation era. 

In Fig. [ll] we compare the evolution of S c for these two models. It is evident that the suppression of S c in the scalr 
field model starts much earlier than in the MDM model. This is also clear when we compare n e ff in the two cases. 

In conclusion, all these effects we have just discussed combine to give additional suppression in the </>CDM model 
as compared to the MDM model. 



E. The temperature anisotropy angular power spectrum 

To complete our analysis of the evolution of perturbations in this cosmology, we shall now look at the effect of the 
scalar field on the cosmic microwave background. The temperature anisotropy measured in a given direction of the 
sky can be expanded in spherical harmonics as 

— (n) = ^2a hn Y tm (n) (53) 

We work again in the framework of inflationary cosmologies so the ai m s are Gaussian random variables which satisfy 

(a}, m ,ae m ) = Ce5g<eS m i m (54) 

The angular power spectrum, Cg , contains all the information about the statistical properties of the cosmic microwave 
background. One can relate it to the temperature brightness function we derived above through 

C t = A-K ^|A Tf (fc,r )| 2 (55) 

where we assume again that we are considering the scale-invariant, inflationary scenario. 
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FIG. 10. The evolution of S at horizon crossing (k 
dark matter cosmology (dashed line). 



5hMpc for the scalar field cosmology (solid line) and for the mixed 




FIG. 11. In the left panel we plot the evolution of 8 C (k = IMpc - ) for the scalar field cosmology (solid line) and the mixed 
dark matter cosmology. In the center panel we plot the suppression factor relative to the standard cold dark matter scenario. 
In the right panel we plot the dimensionless growth rate for these two modes. 
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FIG. 12. The angular power spectrum of CMB anisotropies for different values of These power spectra are not normalised 
to COBE. 



In ref [p2| , the authors presented a useful simplification with which one can understand the various features in the 
angular power spectrum. We shall present a simplified version, in the synchronous gauge and use it to understand 

the effect will have on the CgS. To a reasonable approximation we can write (defining \ ~~ h+6 ^ x 



2fe 2 



A«(J;,to) 



2x}(k,T*)je{k(T -n)) 
I 

ji-i(k(T - r*)) 



1 



L 2i + r 



21+1 



ji-i(k(T - r»))] 



+ / (j)(fe,r) + x(*,T))ji(fe(7TD-r)) 



(56) 



assumin g in stantaneous recombination at t* (we do not consider Silk damping in this discussion). 

In Fig h2 we plot the angular power spectrum of the CMB for a few values of Cg. Comparing to a SCDM Ci we see 
two effects. Firstly, for £ > 200 we see that the peaks are slightly shifted to smaller scales. In [j5l| a clear analysis was 
presented which explains the location of the peaks in the MDM model scenario, and we will now see that the same 
effect is present here, albeit in the opposite direction. The anisotropies at these scales are generated primarily at r* 
and for the purpose of this discussion we can assume that the first term in (|5^) dominates. As shown in J52]] the peak 
structure of this term is given by 



5 7 + 2\ oc cos kr s (r* ) 



(57) 



Now one of the effects described in 



where r s is the sound horizon in the baryon-photon fluid, r s (r) = 3 [ 1+ ^( T /) ] • 
and which we see here is that for a given redshift (during a period of time which contains the time of recombination) , 
r s is different in the presence of the scalar field than in it absence. This is to be expected as deep in the radiation 
era we now have three components contributing, the radiation (whose energy density is set by Tqmb, the massless 
neutrinos (whose energy density is set by T u ) and the scalar field whose energy density is a constant fraction of the 
total energy density. Therefore, for a given redshift, the total energy density in the presence of the scalar field is larger 
than the total energy density in the absence of the scalar field. As a consequence the expansion rate will be larger 
and the conformal horizon will be smaller. In Fig. [l3] we plot the ratio of energy densities and the ratio of comoving 
horizons as a function of scale factor. At recombination (a ~ 10~ 3 ), the conformal horizon (and consequently the 
sound horizon) is smaller, and as seen in Fig. [l2]the peaks are shifted to the right. 
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FIG. 13. The top panel shows the ratio of the total energy density in a universe with 0.$ = 0.1 to that in a scalar field free 
universe. The bottom panel shows the ratio of conformal horizons. 
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FIG. 14. A comparison of the source of acoustic oscillations in <5 7 for k = O^hMpc^ 1 for SCDM (dashed) and 0CDM (solid) 
with 0.4, = 0.1 and H = 50Kms _1 Mpc _1 . 
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FIG. 15. A comparison of the angular power spectra of temperature anisotropies for five COBE normalized models with the 
current experimental situation. The models are sCDM in solid, (h = .5), ACDM in long dash, (Qa = -6, h = .65), MDM in 
dotted, (fi„ = .2 and h = .5), 0CDM in dot-short dash, = .08, h = .5) and 0CDM in dot-long dash, = .12, h = .65). 
All of them have n b h 2 = 0.0125. 



Finally, the other effect f2^ has on the peaks is to boost their amplitude. Let us assume c 2 
the radiation era. As shown in the section III C we can write 



5 k \ 



3 



i.e. we are deep in 



(58) 



so the acoustic oscillations are sourced by S c . One finds that this source is larger in the presence of the scalar field 
than in its absence (see Fig. [l4]) . The increase in the amplitude of this driving term will increase the amplitude in i5 7 
and lead to the increase by a few percent of the acoustic peaks of the angular power spectrum. 



IV. CONSTRAINTS FROM COBE AND LARGE SCALE STRUCTURE 



In the previous section we have discussed in some detail the evolution of density perturbations and temperature 
anisotropies in the </>CDM cosmology. This approximate analysis indicated that results should be similar to those in 
the MDM model, except that there should be additional suppression of power in A 2 (k) on small scales. Using the 
full results of our numerical evolution we now compare our model with the observational constraints which measure 
fluctuations on a wide range of scales. We first compare the Cgs to the COBE data and calculate the normalization 
of the perturbations model for different Q^. For completeness we plot a selection of Ce compared to a compilation of 
data sets and other candidate theories. We then compare A 2 (fc) with the observational A 2 (fc) rendered in |M] from 
a compilation of surveys. For this section we define ^CDMi to be a universe with = 0.08, Ho = 50kms _ MVIpc _1 
and 0CDM2 to be a universe with = 0.12, Ho = 65kms _1 Mpc _1 . 

The past five years has seen a tremendous growth in experimental physics of the CMB. Over twenty experimental 
groups have reported detections of fluctuations in the CMB and a rough picture is emerging of the angular power 
spectrum. It is fair to say that the most uncontroversial and useful measurement that we have is that of COBE, 
which tells us that on scales larger than 10° the fluctuation are approximately scale invariant with a Q rm s = ISfJ-K. 
Measurements on smaller scales seem to indicate a rise in the power spectrum, but a convincing constraint is still 
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FIG. 16. The left panel shows the small I angular power spectrum for a family of </>CDM models. The models are Q$ = 0, 
0.04, 0.08 and 0.12 in order of increasing "tilt". The right panel shows the maximum likelihood of the best fit model as a 
function of f2</,. We have fixed h = 0.5 and Qth 2 — 0.0125. 



lacking. In Fig [lj we present a compilation of measurements of |B3[ as compared to two </>CDM models and a few 
candidate rival models. Clearly there is still is a large spread although an overall shape is emerging. 

In the previous section we described the effect that fi^ would have on the CVs. For small is the dominant effect to 
note is the increase in power of the acoustic peaks relative to the large-scale, scale- invariant plateau. The larger ia Q^,, 
the larger is the boost and therefore the smaller the is which are affected. In practice it introduces an effective "tilt" 
in the large angle power spectrum as can be seen in Fig It is useful to quantify how good a fit (/>CDM Ces are to the 
COBE data. The correct framework to work with is maximum likelihood analysis. In fl54| the authors have supplied 
us with an efficient way of evaluating the likelihood of a given mode relative to purely scale invariant fluctuations. 
On the left panel we plot the likelihood of the best fit model as a function of f2^. There are two important things to 
note. Firstly the well known fact that, if one includes the quadrupole, the COBE data favours more tilted models 
and therefore a larger fi^. Secondly, the likelihood function for should be very flat; indeed the dependence of the 
tilt on fi<£ is very weak. As explained in the previous section the effect of on the Ci s is of order a few percent and 
concentrated at large is. There will be little variation on COBE scales. 

One of the key observational constraints for any class of models is the mass variance per unit interval in In k as 
defined in Eq. |4^. This quantifies the amount of clustering over a range of scales. In |7| the authors compiled a 
series of surveys and attempted to extract what they believe to be the underlying A 2 (fc) of the linear density field. 
This involved a series of corrections: Firstly, the assumption that the different samples were biased in different ways 
with respect to the underlying density field; secondly, that there are redshift distortions in the observed structures; 
and finally, that some of the structures have undergone non-linear collapse. This final correction is model dependent 
and in principal great care should be taken in making definitive comparisons between our theoretical A 2 (fc) and that 
presented in In practice we shall assume that possible corrections are small and compare them. In a future 
publication we shall analyse the non-linear features of this theory. 

In Fig. |j~7| we p lot a family of COBE normalized A 2 (k) with h = 0.5 and Q,bh? — 0.0125. We clearly see the features 
described in IIIB, i.e. the larger the £1^, the smaller the k for which A 2 (k) departs from scale invariance. In the other 
panel we can also see how it differs from a MDM model with the same background cosmological parameters. For the 
same energy density in exotic matter component (i.e. <f> or massive v) there is more suppression in the 0CDM case. 
Finally we see how it compares to the data of 0. We find that, for D,^ in the range 0.08 — 0.12 we can match the 
data with as good agreement as the MDM model and some other candidate models. This is displayed in Fig [l^. In 
the following table we tabulate the \2 values (with 15 degrees of freedom) of these models, in increasing goodness of 
fit. ifiCDM is competetive with the best fit model of MDM. 
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FIG. 17. On the left panel we plot the dependence of A 2 (fe) = -jj^j; on f2^>. For h = .5, Q b h 2 = 0.0125 we show (in order 
of decreasing amplitude) plots for Q,^, — 0, .02, .04, .06, .08, .10 and .12. In the right panel we compare a family of MDM 
models (dashed, with Q, v = 0.04, 0.08 and 0.12 in order of decreasing amplitude) with a family of 0CDM models (dashed, with 
0,$ = 0.04, 0.08 and 0.12 in order of decreasing amplitude). 
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FIG. 18. A comparison of the mass variance per unit interval in lnfc, A 2 (fc) = if , , for five COBE normalized models with 
a rendition of the linear power spectrum from various data sets (corrections for non-linearity, redshift distortions and biasing 
have been introduced ). The models are sCDM in solid (h = .5), ACDM in long dash (£l\ = .6, h = .65), MDM in dotted 
(Q„ = .2 and h = .5), (f>CDM in dot-short dash (Slj, = .08, h = .5) and (f>GDM in dot- long dash (f^ = .12, h = .65). All of 
them have Q b h 2 = 0.0125. 
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One can fit A 2 (fc) for these models (to 10% in the range < < 0.16) with: 



A 2 (fc) =D(k,n^)A 2 (k) CDM 
= (1 + 1.5^-10^) 



1 + 5fc 11295 + 4fc 2 
1 + 1 x Wk 2 - 259 



2. 259 x 1.1501' 



(59) 



where A 2 (k)cDM is the COBE normalized CDM mass variance, and contains all the dependence on remaining 
cosmological parameters such as Hq and fib |55| . 

A useful quantity to work with is that characterising the mass fluctuations on i?h _1 Mpc scales 



a 2 (R) 



^A 2 (fc)( 3jl(fcJ?) ) 2 . 
k kR 



(60) 



In particular it has become the norm to compare this quantity at 8h _1 Mpc with the abundances of rich clusters. It 
is premature to use this constraint with the our current understanding of 0CDM. The best measurements of such 
abundances involve an estimate of the number density of X ray clusters of a given surface temperature. To relate 
these temperatures to masses in an accurate way one has to rely on N-body simulations of clusters. This has been 
done for a few cosmologies and we can use the results they use as a rough guide but care should be taken with using 
such results at face value. The fact that we can fit A 2 (fc) to the data of M is already strong indication that we are 
on the right track. If we use the values of |5l| we have erg ~ 0.5 — 0.8. A good fit to a$ is 



cts(^) 



-8.7Qi 



CDM 



(61) 



where <x£ dm is the COBE normalized sCDM erg . As would be expected, for the range of values for which we get a 
good agreement with |?J, we also match the cluster abundance constraints. 

Finally it is desirable to make a comparison with some measure of small scale clustering at early times. In fl59[ , the 
authors used a simple analytic estimate of the fraction of collapsed objects at redshift z — 3 and z — 4 to show that for 
f2i, = 0.3, MDM models predict too little structure as compared to that inferred from the Lyman-a measurements JS0| . 
More recently in |61|, the authors considered a larger range of cosmological parameters and found that constraints 
from the Lyman-a systems could be sufficiently restrictive to rule out a large range of models. From Fig [l8] we can 
see that 0CDM should fare better than MDM on very small scales. This is easy to understand: We argued in Section 
III that effectiveness of (/>CDM was mainly due to the fact that the scalar field free-streaming scale grows with time 
while the massive neutrino free streaming scale decays with time. We then need a larger amount of massive neutrinos 
to fit both COBE and the cluster abundances in the MDM model than the amount of scalar field in 0CDM. On much 
smaller scales (the scales probed by Lyman-a systems) i.e. scales smaller than the massive neutrino free streaming 
scale, perturbations in MDM should be more supressed than in </)CDM. This means 0CDM should fare better than 
MDM with regards to the Lyman-a constraints. A preliminary check on our model can be done using the technique of 
[jsll . In brief one can make a conservative estimate that a fraction f gas of matter is in gas at that time and set bounds 



on the amount of objects with masses greater than 10 10 (1 
can then derive a bound: 



n (j ,)- 1 h- 1 M Q . Using the Press-Schecter formalism one 



erfc( 
erfc( 



1.7 



V2a(R, z) 
1.7 

V2a(R,z) 



)>0.16h/f gas z = 3 



) > O.lOAh/fg 



z = 4 



(62) 



where erfc is the complementary error function and R = 0.1 — 0.2Mpc. For 0CDM in the range of h and that 



we have been considering we find that if f g 



1, these models are consistent with this constraint. Note that the 



MDM model already has serious problems with this conservative constraint. If we consider a less a conservative 
constraint and take f gas ~ 0.1 as seems to be indicated by hydrodynamical studies, than 0CDM is inconsistent with 
these measurements. More detailed observations and modeling of Lyman-a systems will supply us with a very strong 
constraint on this class of models. 
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FIG. 19. The shaded region represents the allowed region of parameter space consistent with the Lyman-a constraints of ]60| . 
The solid hatched region corresponds to the values of h and $7^ excluded for f gas = 1- The dotted hatched region corresponds 
to the values of h and Q$ excluded for f gas — 0.1. 



V. CONCLUSIONS 



In the first part of this paper we gave the motivation for considering the particular scalar field cosmology we have 
now studied in detail: The addition made to the standard SCDM model does not involve any tuning of the type 
involved in other modifications, in that no energy scale characteristic of the universe at recent epochs is invoked. 
The form of the required potential is one which arises in many particle physics models, with values of the single free 
parameter of the order required. Having analysed the model and determined the best fit to structure formation, we 
now first co mmen t on these aspects of the model. 

In Section HE we argued that, in typical inflationary models with the usual mechanism of reheating, one would 
expect the attractor solution to be established well prior to nucleosynthesis. In this case we should therefore compare 
our best fit to the constraint from the latter, as given in (p0|), which is satisfied (clearly for the conservative bound, 
and marginally for the tighter one). Further, in this case we need to check that the value Q,bh? — 0.0125 which is 
equivalent to the baryon to photon ratio 7710 = 3.3 (in units of 10~ 10 ) is within the allowed range at the corresponding 
expansion rate. This is not a question which is easy to answer simply because most nucleosynthesis calculations give 
results only in terms of the range of 7710 allowed at the standard model expansion rate, or the maximum additonal 
number degree of freedom allowed. For the conservative criteria ]4l]j43] | which gave the weaker upper bound in ( pp| ) 
the allowed range narrows from rj W £ [1.65, 8.9] at = to the lower bound 7710 = 1.65 at fL, = 0.15. Extrapolating 
the use of these criteria for a more restrictive case for which the required data is given in |l2| (i.e. the allowed range 
of rjiQ is given as a function of the energy density in an extra component) the value 7710 = 3.3 would appear to be in 
the allowed range at = 0.1. For a more restrictive set of nucleosynthesis constraints a slightly lower value of ilbh 2 
might be required for consistency. 

In an alternative model of reheating (which is by construction associated with the existence of this same type of 
potential), we saw that the time of re-entry to the attractor could be after nucleosynthesis, at a time which depends 
on both and the other parameter in this model Hi, the expansion rate at the end of inflation. By making 
the assumption that the attractor is established prior to nucleosynthesis we restricted ourselves to a (large) part of 
parameter space. From Fig || we can read off that this corresponds, for f2^ w 0.1, to Hi > 10 14 GeV or p^ 4 > 10 16 
GeV. As we noted at the end of section this is consistent with what would guess would be the most natural range of 
these parameters in this model. On the other hand, it is inconsistent with models in which the phase of scalar field 
domination continues until just before nucleosynthesis (with the consequences described in [ flB") ) since they correspond 
to the line defining the lower bound from nucleosynthesis in Fig ^. For entry to the attractor before today, it can be 
see from Fig that one requires in this case Q$ > 0.22. It would be interesting to analyse the effect on structure 
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formation in these models, and indeed in all of the parameter space for this model excluded by us in our present 
analysis. 

In terms of A our best-fit corresponds to the range A £ [5, 6.1]. As we discussed in Section [IB this value is certainly 
of the order observed in the fundamental particle physical theories of which they have been observed to be a generic 
feature, and may even be in the precise range found in certain theories. This suggests the exciting possibility of 
ultimately linking the cosmological features which would provide a signature for these fields to details of physics 
at the Planck scale. From the point of view of particle physics motivated model-building it would be particularly 
interesting also to look at models where the simple exponential potential represents the asymptotic behaviour of a 
potential which can support inflation in another region, since this would be likely to produce a very constrained model 
(with reheating as discussed). 

In Section 111 we analysed in detail the evolution of perturbations in 0CDM. We found that perturbations in the 
scalar field on subhorizon scales decayed, leading to a suppression of power on small scales. The similarities with 
MDM led us to pursue the comparison in more detail. We found that the contribution from the scalar field was more 
efficient at suppressing perturbations in the CDM than massive neutrinos. We showed that this was due to a simple 
difference in the evolution of the "free-streaming" scale in the two theories: In c/>CDM the free streaming scale grows 
with the horizon while in MDM the freestreaming scale decays as 1/r. This means that perturbations in the scalar 
field never grow once they come into the horizon, in contrast to perturbations in the massive neutrino which end up 
clumping after some finite time. We analysed the effect the scalar field would have on the CMB and found it to be 
small but distinctly different from that of MDM. 

In Section [fy| we used the results of a Boltzmann-Einstein solver to test how well this class of models fared when 
compared to various astrophysical data. Because of the weak effect the scalar field has on the CMB, the angular power 
spectrum is effectively (with the current accuracy of experiments) as in SCDM. Using the COBE data we normalized 
these theories and compared the mass variance per logarithmic interval in k to the one estimated in |7j]. Our models 
fared as well, or better, than competing flat universe models. A comparison with an estimate of the mass variance at 
8h _1 Mpc from the abundances of rich clusters gave the same results. We finally compared the amount of structure 
at high redshift our model predicts, as compared to that infered from Lyman-a systems. This has proven to be a 
serious problem for MDM models. We saw that our model is consistent, albeit marginally, with these constraints. 

Lastly a few further commments on other related issues which it would be interesting to investigate: 

(i) We assumed an initial flat adiabatic spectrum of perturbations, in line with the most generic type of inflation. 
Within the context of inflation one can of course have different spectra etc., and within the context of some very well 
motivated form for the inflationary part of the potential in the alternative reheating model, it would be interesting 
to look at the combined effect on structure formation. In more general, an interesting feature of the exponential 
which would be worth investigating is the fact that the attractor is also an attractor for isocurvature fluctuations, and 
hence the assumption of adiabatic initial conditions might have a much more general motivation than the standard 
inflationary one. 

(ii) Inflation was assumed simply because it is the paradigmatic model. The existence of the exponential scalar 
field might of course have an effect on any cosmology. In an open cosmology, for example, the analagous attractor 
also exists in the curvature dominated regime and the asymptotic state has a scalar field energy scaling as I /a 2 . 

(iii) We have shown that the effect of the exponential scalar field on the angular power spectrum of the CMB is 
quite small, much like the case of MDM. However, recent high precision analysis of parameter estimation from the 
CMB (as one would expect from the satellite missions) indicates that one may achieve a precision of or 
better. This opens up the interesting possibility of actually trying to constrain ft^ with the CMB. 

(iv) Besides the issue of consistency with entry into the attractor prior to nucleosynthesis which would motivate 
the study of the dependence on fl},h 2 away from the SCDM value we assumed, there are further observational reasons 
for doing so. The recent measurements of fl62| indicate that one may have a higher baryon content than previously 
expected. Since a higher baryon content leads to less structure on smaller scales, a best fit to large scale stucture 
constraints would be obtained with a smaller value of fi^. The effect on the Lyman-a constraints might be more 
dramatic, since there would be a competition between the suppression of power and the increase in the amount of gas 
simply due to the fact that there are more baryons around. That the result is not immediately evident can be seen 
from the analysis of MDM in |)si| . Further, to determine whether entry to the attractor prior to nucleosynthesis is 
consistent in this case, one would have to determine whether this decrease in $7^ would broaden the allowed range in 
r) sufficiently. 

(v) One of the main problems for the MDM model is the overwhelming evidence for structure at high redshift. 
Critical universe models with massive neutrinos typically underproduce structure as probed by Lyman-a systems 
and high redshift cluster abundances pty . This problem is somewhat alleviated in the case of (/>CDM: The different 
evolution of the free-streaming scale in </>CDM leads to more power on small scales (and consequently high redshifts) 
relative to MDM. This certainly provides strong motivation for further study of the evolution of perturbations on 
small scales and high redshift in these models. 
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From the point of view of structure formation we have described a new model which has the same qualitative 
features as MDM. Unlike other scalar field cosmologies, which affect the local expansion rate (up to redshift of a 
few), our background evolution is exactly matter dominated, the modification arising at the perturbation level. Given 
the wealth of current data on small scales, at recent redshifts, this raises the question of whether this is the right 
approach to the construction models of structure formation: The existence of structure at high redshift combined 
with the small scale velocity dispersion today at lh _1 Mpc seems to argue for a strong modification of the growth 
rate of perturbations in the last few redshifts. This would point towards a low density universe. However, until we 
have a more detailed understanding of the non-linear evolution of perturbations in models such as 0CDM (on scales 
between 0.1 and 10h _1 Mpc), these models should not be ruled out. We are currently analysing the non-linear regime 
of 0CDM using an N-body code. 
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